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ADMINISTRATIVE  INFORMATION 

This  program  was  authorized  end  funded  by  the  Nava)  Ship  Engineering 
Center  Wbrk  Request  Number  N65197-75-P0-5131*  and  is  identified  as  DTNSRDC 
Work  Unit  Number  1-1560-012. 


INTRODUCTION 

A ship's  response  in  irregular  waves  may  be  derived  from  the  irregular- 

wave  spectrum  and  the  regular-wave  response  by  the  linear  superposition 

theory  of  St.  Denis  and  Pierson.*  The  regular-wave  transfer  function 

may  be  obtained  from  experimental  data  or  from  analytical  predictions  provided 

2 3 

by  existing  computer  routines.  In  the  past,  idealized  wave  spectra  such 
as  the  Pierson-Hoskowitz  formula  for  fully  developed  seas  and  the  Bretschneider 
two-parameter  sea  formula  have  been  used  to  represent  irregular  waves  of  various 
significant  wave  heights. 

In  recent  years,  DTNSRDC  has  obtained  the  data  for  323  point  spectra 

with  significant  wave  height  from  0.*i  to  kk.S  feet  which  were  measured  at 

4 

Station  INDIA  in  the  North  Atlantic.  More  recently,  a series  of  mathematical 
wave  spectra  involving  six  parameters  has  been  developed  by  Ochi'*  in  an 
attempt  to  systematically  represent  the  variety  of  wave  spectra  shapes  which 
may  occur  for  any  given  significant  wave  height- 
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St.  Denis,  M.  and  W.J.  Pierson,  "On  the  Motions  of  Ships  in  Confused  Seas," 
SNAHE  Trans.  Vol.  62  (1951*) 

Frank,  W.  and  N.  Salvesen,  ,rThe  Frank  Close-Fit  Ship-Motion  Computer  Program," 
NSRDC  Report  3289  (1970) 

^Meyers,  W.G.,  D.J.  Sheridan,  and  N.  Salvesen,  "Manual,  NSRDC  Ship-Motion  and 
Sea  Load  Computer  Program,  "NSRDC  Report  3376  (1975) 

4 

Miles,  M. , "Wave  Spectra  Estimated  from  a Stratified  Sample  of  323  North 
Atlantic  Wave  Records,"  N.R.C.  Report  LTR-SH-118-A  (May  1972). 

**0chi,  M.K.  and  E.N.  Hubble,  "On  Six-Parameter  Wave  Spectra,"  Proceedings 
of  the  Fifteenth  Conference  on  Coastal  Engineering  (July  1976)- 


This  report  describes  a computer  program  for  deriving  long-crested 
irregular-wave  responses  by  the  method  of  linear  superposition  using  the 
Station  INDIA  sea  spectra  and/or  the  six-parameter  wave  spectra.  The 
Pierson-Hoskowitz  spectra  may  also  be  used  for  comparison.  The  program, 
entitled  SRIW,  is  written  in  FORTRAN  for  operation  on  the  CDC  6700 
computer  at  OTNSRDC.  Cal  comp  plots  of  the  results  may  also  be  obtained. 

PROCEDURE 

Any  irregular-wave  response  — pitch,  heave,  roll,  yaw,  sway,  surge, 
relative  motion,  forces,  moments,  etc-  — may  be  calculated  if  the 
regular-wave  transfer  function  or  response  amplitude  operator  is  available 
in  Gi?e  c-f  the  formats  specified  in  the  Appendix.  The  program  computes 
significant  values  of  the  response,  most  probable  extreme  values,  and 
extreme  values  for  design  consideration.  In  addition,  responses  such  as 
relative  velocity,  deck  wetnesses,  slams,  etc.,  may  be  derived  from 
relative  vertical  motion.  Vertical  and  lateral  accelerations  may  also 
be  derived  from  absolute  vertical  motion  and  roll,  respectively.  Results 
are  tabulated,  and  plotted  if  desired,  cs  a function  of  significant  wave 
height.  A sample  plot  is  shewn  in  Figure  1. 

SIGNIFICANT  VALUES 

Variances  and  significant  values,  i.e.,  average  of  highest  1/3*  are 
derived  in  the  following  manner.  Integration  is  performed  in  the  wave 
frequency  u domain  with  limits  of  0 to  2.0  rad/sec  as  suggested  in 
Reference  6. 

Hubble,  E-h.,  and  J.B.  Hadler,  "Prediction  of  Ship  Motion  in  Regular  and 
Irregular  Head  Seas,"  NSRCC  Report  SPD-623_01  (April  1975) 
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(1)  Wave  frequency  in  rad/sec:  u ■ /Zitg/X 

2 

(2)  Encounter  frequency  in  rad/sec:  uft  • u - u cos  y V/g 

2 

tere  g « gravitational  acceleration  in  ft/sec 
X * wave  length  in  ft 
V * ship  speed  in  ft/sec 

y * heading  angle,  180  deg  for  head  seas,  90  deg  for  beam  seas 
0 deg  for  following  seas 


-1  + J\~- 


u = 


4 cos  y — ue 

for  90°  < y < 270° 


-2  cos  y — 
9 


(3)  Response  amplitude  operator  at  V,  u . and  y transformed  to  u: 

e 2 

R A 0 * x (u)  = (rc9ular-%*avfe  response  amplitude. 

' wave  amplitude  in  ft  ’ 


(4)  irregular-wave  spectral  ordinate  at  u in  — : 

rad/sec 

S (to) — Values  from  Reference  4 for  Station  INDIA  Spectra 
S(u) — Values  derived  from  equation  and  constants  in  Reference  5 
for  Ochi  six-parameter  spectra 

(5)  Variance  of  displacement: 


= n?Q  = J2m0  x(u 


X (u)  S (u)  du 


(6)  Variance  of  velocity: 
•2.0 


E = m 
v 2 


-r°  j 
■Jn 


X(u)  S (to)  du 


(7)  Variance  of  acceleration: 

2.0 


rl.O 

‘"‘•'Jo  “« 


X(u)  S(u)  du 


(8)  Significant  value  of  the  response,  i.e.,  displacement: 
X,  = 2.0025  »C 


Pi|SSi 


r-  3HH*,*sP^feRV; 


(9)  Significant  velocity:  V^  “ 2-0025 

(10)  Significant  acceleration:  Aj^  “ 2-0025  /E^ 

EXTREME  VALUES 

Extreme  values  are  derived  from  the  following  relationship  given  by  Ochi 

(11)  Extreme  value  of  the  response: 


: 7 


where  a = 1.0  for  most  probable  extreme  value 

a = 0.01  for  extreme  value  for  design  use 
with  99  percent  assurance 
T = maximum  duration  time  in  hours  for  wav°s  of 
given  wave  height  (see  Figure  2,  extracted 
from  Reference  8) 

SUPPLEMENTARY  RESPONSES 

Additional  statistical  responses  may  be  derived  from  relative  vertical 

g 

motion  R.A.O.'s  as  follows: 

(12)  Number  of  deck  wetnesses,  water  contacts  on  catamaran 
cross-structure,  or  hydrofoil  broachings  per  hour: 


N - 3600 
c ~ 2ir 


& 


d eV/(2Ed) 


70chi,  M.K-,  "On  Prediction  of  Extreme  Values,"  Journal  of  Ship  Research,  Vol.  17, 
No.  I (1973) 

®0chi,  M.K. , and  Hotter,  L.E.,  "Prediction  of  Extreme  Ship  Responses  in  Rough 
Seas  of  the  North  Atlantic,"  Proceedings  of  the  Symposium  on  the  Response  of 
Marine  Vehicles  and  Structures  in  Waves,  Univ.  College  London,  (197*») 

^Ochi,  M.K. , "Prediction  of  Occurrence  and  Severity  of  Ship  Slamming  at  Sea," 
Proceedings  of  the  5th  Symposium  on  Naval  Hydrodynamics,  Office  of  Naval  Research 
(196A)  k 
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c. 


^ JV* 
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ei  ■ 

r* 


ft 


r<v4^Sry,-i>^^ 


where 


Ej  » variance  of  relative  vertical  displacement  in  ft 


Ey  ” variance  of  relative  vertical  velocity  in  (ft/sec)i 


C * absolute  vertical  distance  in  ft  from  water  surface 
o 


(running  in  still  water)  to  appropriate  structure 
(13)  Humber  of  slams  per  hour 

N . 3|»0  X7T  e-(C02/(2Ed)  - f42/(2Eu)] 
s 2tt  I v d 


where 


r-j.  = threshold  relative  velocity  in  ft/sec 


( 1 4*)  Significant  pressure  on  catamaran  cross-structure  in  psi: 
P, 


1/3  =2.1k  (2Ey) 


where 


k = coefficient  dependent  upon  cross-structure  section 
shape 


PERCENTAGE  EXCEEDANCE 

Cumulative  statistics  are  applied  to  the  Station  INDIA  data,  if  desired, 
to  predict  the  probability  of  occurrence  of  the  response  over  a one-year 
period  in  the  North  Atlantic.  Since  the  Station  INDIA  spectra  do  not 
necessarily  represent  an  average  of  the  real  ocean  environment,  each  spectra 


is  given  a weighting  factor.  These  weights  were  derived  from  meteorological 

10 


data  presented  by  Hogben  and  Lumb  for  Area  2 which  represents  the  Eastern 
half  of  the  North  Atlantic  and  encompasses  Station  INDIA.  The  results  are 
presented  as  percentage  exceedance  of  the  response.  Sample  plots  are  shown 
in  Figure  3- 


10, 


Hogben,  N.  and  F.E.  Lumb,  "Ocean  Wave  Statistics,"  (19&7) 


5 


RELATIVE  BOW  MOTION  IN  FEET 


SIGNIFICANT  VALUES 


HEADING=180 -0  DEG 


L=220  *0  FT 


V=1U.0  KNOTS 


o 
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Figure  1 — 


Sample  Plot  of  Irregular-Wave  Response  versus  Significant  Wave  Height 
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CONTROL  CARDS  FOR  SRTW  AT  NAVSEC.  USING  ON-LINE  1700  CALCQMP  PLOTTER 


XXXXSRI.Cm62000«T200«P3. 

CHARGE  . X X XX  • J.l  J J J J J J J J » PR  • 

HAP. OFF. 

ATTACH . SRI W • SH I ORESRONSE I RREGUL ARHA VES . I D.CHNH . 

ATTACH.TAOE2.IN0lASPECTRA.I0.CHNH, 

ATTACH. TAOE3.0CHISPECTRA.ID=CHNH, 

ATTACH  . T AOE^»  • PMSE  AFORMUL  A • I D*CHNH  . 

ATTACH* CALClTO. 

attach.calcfn, 
ldset.li a=C ALC 1 70 . 
ldset.li r=calcfn. 

SRIW. 

ROUTE • PLOT • OC=PT . TI D= AF . 

oooooooooooooooooooooo  end  of  record 

INPtlT  DATA  CARDS  ' 

oooooooooooooooooooooo  end  of  FlLF 


1 


J 

‘I 

1 

:J 


CONTROL  CAROS  FOR  SRlW  AT  uTnSRDC.  USING  CALCQHP  936  PLOTTER 


XXXXSPI »CM<>2000.Ti?00.MTl.P3. 

CHARGE  * XX  XX . J J J J J J J J JJ.PR. 

MAP. OFF. 

VSN.TAPF1=CAOOOO. 

RE0UEST.TAPE7.HI.RING. 

ATTACH . SRIW* SHI PRESPONSE iRRtGULARHA VES . I D*CHNh . 
ATTACH. TAPE2. IN3IASPECTRA. IU=CHNH. 

ATTACH. TAPE3.nCHlSPLCTRA,ID=CHNH. 

ATTACH. TAPE* • PMSE AFORMUL A , I D=CHNH. 

ATTACH. CALCR3*. 

ATTACH.CALCFN. 

LOSET.LIB=CAlC<}36. 

LOSET,LIR=CALrFN. 

SRIK. 

OOOOOOOOOOOOOOOOOOOOOO  END  OF  RECORD 

INPUT  DATA  CARDS 

OOOOOOOOOOOOOOOOOOOOOO  end  of  file 


CONTROL  CARDS  FOR  SRTW  WITHUUT  ANY  PLOTS 
XXXXSRl .CM62000.T100"P3. 

charge . xxxx . jjjjjjjj jj.pr. 

map. OFF. 

ATTACH.SRIH.SHlPRESPONSEIRRfcGULARHAVES.IO.CHNH. 

ATTACH.TAPE2»TNDIASPECTRA,in=CHNH. 

ATTACH. TAPE 3 .OCHl SPECTRA. ID=CHNH. 

ATTACH  . T APE-*  . PMSE  AFORMUL  A . I D=CHNH  . 

SPIW. 

OOOOOOOOOOOOOOOOOOOOOO  end  of  record 

INPtlT  DATA  CAROS 

OOOOOOOOOOOOOOOOOOOOOO  end  of  file 
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CARD 

COLUMNS 

FORMAT 

FORTRAN 

EXPLANATION 

1 

1 - 3 

13 

NSHIPS 

Number  of 

cases  to  be  run;  no  limit 

1 

11-40 

3A10 

NAME 

User's  nan 

e,  code,  phone  number 

2 

1 -4o 

4A10 

SHIPN 

Main  title 

:,  such  as  ship  name,  etc. 

2 

41-70 

3A10 

RESP 

Label  for 

irregular-wave  response  to 

conputed,  Including  dimensions. 

Examples: 

PITCH  AMPLITUDE  IN  DEGREES 
R.B.M.  AT  STATION  1.5  IN  FEET 
Set  IRM  to  1 if  the  response  to  be  conputed 
is  relar.ve  notion.  Otherwise,  leave  blank 
or  set  IRM  to  0. 


1-3  13 


IFORM  Code  to  identify  form  in  which  the  transfer 
function  is  input 


IF0RM=1,  if  TF  = (i){£)  and  (|)2  = (TF*C)Z 
IF0RM=2,  if  TF  = (f ) (~~)  and  (|)2  =(TF=C*^) 


I F0RM=3 , if  TF  = (i)(|)2  and  (|)2  =(TF*C) 


g ra 

^ "L  S 


Definitions: 


TF  * input  transfer  function 

X - response  amplitude  in  same  dimensions  specified  on  Card  2 
R * regular-wave  amplitude  in  feet 

X = regular-wave  length  in  feet 

K * wave  number  * 2n  /X 
KR  = wave  slope  in  radians  * 2 if  R/X 

C = constant,  input  on  Card  k 
X 2 

(— ) = response  amplitude  operator 


EXAMPLES  OF  POSSIBLE  TRANSFER  FUNCTION  FORMS 


Irregular-Wave  Response  input  Transfer  Function  (TF) 


Heave  in  feet 


Heave  amplitude  in  feet 
Wave  amplitude  in  feet 


Moment  in  ft-tons 
Pitch  in  radians 
Roll  in  degrees 
Yaw  in  degrees 


Moment  amplitude  in  ft-lbs  x 10  ^ 
Wave  amplitude  in  feet 

Pitch  amplitude  in  radians 
Wave  slope  in  radians 

Roll  amplitude  in  radians 
Wave  slope  in  radians 

2 

,Yaw  amplitude  in  deg, 

'Wave  amplitude  in  ft' 

/Force  amplitude  in  lbs,^ 

'wave  amplitude  in  ft 


I FORM  C Value 
1 1-0 

I ,()6 

1 2250 

2 1-0 

9 360 

z 2n 

3 1-0 


( : ) 

'2240  * 


Force  in  tons 


*vu’ 


tT.  - ^ *”■  - .$3^  'S3*s-*  «•  *»•*  "“"*«.  *>***“■  -ic 


CARD  COLUMNS  FORMAT  FORTRAN 


EXPLANATION 


4-6 


IFREQ  Code  to  identify  frequencies  at  which 


transfer  function  is  input 

I FREQ* 1,  if  exactly  80  values  are  input  corresponding 

to  wave  frequencies  ca  = 0.25,  -05,  -075,  

1-975,  2-0  rad/sec 

IFREQ=2  if  50  valued  or  less  are  input  corresponding 
to  selected  values  of  wave  length/ship  length 
(X/L)  input  on  Cards  6a,  6b,  etc. 

IFREQ=3,  if  50  values  or  less  are  input  corresponding 

to  even  increments  of  nondiraensional  encounter 
frequency  (ca^  l/t/g)  specified  in  Card  6. 

IFREQ=A  if  50  values  or  less  are  input  corresponding 


* ^ 


7-9 


to  selected  values  of  wave  frequency  ta  input 
on  Cards  6a,  6b,  etc. 

Number  of  X/L  values  input  on  Cards  6a,  6b,  etc. 


& 

| 

£ 


10-12 


when  1FREQ=2,  or  number  of  ia  values  input  on 
Cards  6a,  6b,  etc-  when  IFREQ=4.  Maximum 
value  of  NWL  is  50.  NUL  not  used  when  IFREQ=1 
or  IFREQ=2- 

Code  for  supplementing  response  to  be  computed 
in  addition  to  the  basic  response  identified 


on  Card  2. 


ISR  = 0 if  no  supplementary  response  is  required. 


1 


"The  option  IFREQ=3  can  be  used  only  if  90°  < p < 270°  since  this  program  is 
set  up  to  convert  from  ca  to  (a  only  in  this  regime.  See  Equation  (2)  page  3* 
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5. 

> 

SUPPLEMENTARY  RESPONSE 

BASIC  RESPONSE 

VALUE  OF  CLR 

[>  • 

jSR 

TO  BE  COMPUTED 

IDENTIFIED  ON  CARD  2 

REQUIRED 

u 

l 

p 

1 

Significant  Relative 

Relative  Motion  (R.M.) 

- 

J5' 

Velocity  in  ft/sec 

in  feet 

i 

g^ 

2 

Significant  Vertical 

Absolute  Motion  (A.M.) 

- 

i 

Acceleration  / g 

in  feet 

>■ 

jC 

r 

3 

Significant  Lateral 
Acceleration  / g 

Roll  Angle  in  degrees 

CG  to  Deck 

1 

4 

Deck  Wetnesses  per  hour 

R.M.  in  feet. 

WL  to  Deck 

1 

1 

generally  at  the  F.P. 

t 

5 

Hull  Slams  per  hour 

R.M.  in  feet. 

WL  to  Keel 

h 

generally  at  Sta.  3 

fe' 

5 

Water  Contacts  per  hour 

R.M.  in  feet. 

WL  to  C-S 

sr  - 

on  Catamaran  Cross- 

at  fwd  end  of  C-S 

| 

Structure  (C-S) 

* 

7 

C-S  Slams  per  hour 

R.M.  in  ft,  at  C-S 

WL  to  C-S 

j& 

if 

8 

Significant  Pressure 

R.M.  in  ft. 

- 

1 

ji  ■ 

on  C-S  in  psi 

at  fwd  end  of  C-S 

9 

Hydrofoil  Broachings 

R.M.  in  ft. 

WL  to  top 

1 

5S 

per  hour 

at  hydrofoil 

of  foil 

14 


CARD  COLUMNS  FORMAT  FORTRAN 
3 13-15  13  I SEA 

ISEA-1, 
I SEA-2, 
ISEA-3, 

I SEA-4, 

ISEA-5, 

3 16-18  13  IPE 


EXPLANATION 

Code  for  irregular-wave  spectra  to  be  used 

for  Station  INDIA  spectra  only 

for  Ochi  six-parameter  spectra  only 

for  Station  INDIA  spectra  and 

Pierson-Moskowitz  spectra 

for  Ochi  six-parameter  spectra  and 

Pierson-Moskowitz  spectra 

for  all  three  sets  of  spectra 

Set  IPE  to  1 if  cumulative  statistics  are  to 
be  applied  to  the  Station  INDIA  data  to  predict 
the  probability  of  occurrence  of  the  response 
over  a one-year  period  in  the  North  Atlantic. 
Results  will  be  presented  as  percentage 
exceedance. 

If  not  required,  set  IPE  to  0. 


19-21  13 


1PL0T  Code  for  response (s)  to  be  plotted 

1PL0T-0,  if  no  plots  are  required 
IPL0T—1,  if  only  significant  values  of  the  basic 
response  are  to  be  plotted 
1PL0T=2,  if  significant  and  most  probable  extreme 
values  are  to  be  plotted 

!PL0T=3,  if  only  extreme  values  for  design  use  are  to 
be  plotted 

1PL0T-4,  if  only  the  supplementary  response 
designated  by  ISR  is  to  be  plotted 
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EXPLANATION 


CARD 

COLUMNS 

FORMAT 

FORTRAN 

3 

19-21 

13 

1 PLOT-5, 

3 

22-24 

13 

ISYMA 

25-27 


1PLOT-5,  if  all  computed  responses  are  to  be  plotted 


ISYMA  Code  for  symbol  to  be  used  when  plotting 

data  derived  from  INDIA  spectra 
ISYMA-3,  will  produce  the  symbol  + 

There  will  be  one  data  point  for 
each  of  the  323  INDIA  spectra 

ISYMB  Code  for  symbol  to  be  used  when  plotting 

data  for  Ochi  six-paraxeter  spectra 
ISYMB-O,  will  produce  the  symbol  0 

There  will  be  11  data  points  for  each 
wave  height  used. 


3 28-39  13  ISYMC  Code  for  additional  symbol  to  be  used  when 

plotting  the  most  probable  six-parameter  data; 
one  data  point  for  each  wave  height 
ISYMC-11,  will  produce  the  symbol  # 

Since  the  most  probable  six-parameter  data 
is  plotted  twice,  the  resultant  symbol  is  9 
if  I SYMC=0  and  ISYMC=11 

Notes: 

(1)  Responses  are  plotted  as  a function  of  significant  wave  height.  In  addition, 
percentage  exceedance  is  plotted  as  a function  of  the  basic  response,  if  IPE=1. 

(2)  Other  symbols  may  be  used  if  desired;  see  Calcomp  manual. 

(3)  Pierson-Hoskowitz  data  is  plotted  as  a solid  line. 
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CARD 

COLUMNS 

FORMAT 

FORTRAN 

EXPLANATION 

4 

1-8 

F8.2 

BPL 

Ship  length  (Lgp) 

in  feet 

4 

9-16 

F8.2 

DISPL 

Ship  displacement 

in  long  tons  (optional) 

4 

17-2 

F8.2 

VKT 

Ship  speed  (V)  in 

knots 

4 

25*32 

F8.2 

BETA 

Heading  angle  (u) 

in  degrees 

33-4o 


i»l-J«8 


F8.2 


F8.2 


^9-56  F8.2 


57*64  F8.2 


BETA~180,for  head  seas 

BETA»90,  for  beam  seas 

BETA*0,  for  following  seas 

C Constant  for  converting  input  transfer 

function  to  proper  dimensions;  see  IFORM, 
Card  3- 

CLR  Vertical  distance  (absolute  value)  in  feet 

required  for  calculation  of  supplementary 
response;  see  I SR,  Card  3 Distances 
should  be  measured  from  the  running  water- 
line if  possible 

ALPHA  Constant  (a)  required  for  calculation  of 
extreme  values  for  design  use 

ALPHAf.  01 , for  59  percent  assurance 

XKK  Constant  (k)  required  for  calculation  of 
significant  pressure  on  catamaran  cross- 
structure, when  I SR  = 8. 
k is  dependent  upon  the  cross-structure 
section  shape 
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COLUMNS  FORMAT  FORTRAN 


EXPLANATION 


65-72  F8.2  ROOT  Threshold  velocity  (rj  in  ft/sec 

required  for  calculation  of  slams,  when  ISR=5  or 
I SR»7-  If  r±  is  not  input,  program  will 

, • „ /l__  / 520 

assign  a value  of  rx  * 12.0  / BP 


5 1-8  F8.2  YINC(l)  Scale  for  one  inch  on  plots  of  signifi- 

cant values  of  basic  response 

5 9-16  F8-2  Y I NC (2)  Scale  for  one  inch  on  plots  of  most 

probable  extreme  values 

5 17-24  F8.2  Y 1 NC (3)  Scale  for  one  inch  on  plots  of  extreme 

values  for  design  use 

5 25  32  F8.2  YINC (4)  Scale  for  one  inch  on  plots  of  supplementary 

response,  I SR 

Notes:  (1)  If  no  plots  desired,  or  if  program  is  to  set  appropriate  scales, 
leave  Card  5 blank. 

(2)  Abscissa  for  basic  plots  is  significant  wave  height;  scale  for 
one  inch  is  8 ft;  abscissa  is  7 inches  long. 

(3)  Ordinate  for  basic  plots  is  the  appropriate  response  in  dimensions 
specified  on  Card  2,  or  the  fixed  dimensions  specified  for 
supplementary  responses;  ordinate  starts  at  0-0  and  is  8 inches  long. 

(4)  Scales  for  percentage  exceedance  diagrams  are  set  by  program. 
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TRANSFER  FUNCTION  - OPTION  1 


COLUMNS  FORTRAN 


1-8  F8.3  TF(72)  Transfer  function  at 


73~80  F8.3  TF(8l)  Transfer  function  at  w 
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TRANSFER  FUNCTION  - OPTION  2 (IFREQ-2) 


I 


63* 


j.rl 

■sjm 
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CARD 

COLUMNS 

FORMAT 

FORTRAN 

6a 

1-8 

F8.3 

NAVEL (1) 

9-16 

F8.3 

WAVEL (2) 

73-80  F8-3  WAVEL(IO) 


Sb  1-8  F8.3 

Continue  with  10 
values  per  card. 
Max.  is  5 cards. 


EXPLANATION 


»oiwj  ui  — r. i 77 — r 

ship  length  L 

at  which  transfer  function 
is  input 


WAVEL(II)  A/L  can  be  in  either  ascending  or 

• descending  order,  and  need  not 

• be  equally  spaced. 


E 

I 

| 


B 


n 


WAVEL(NVL) 


I 

I 


m 


js* 

I 

►/ 

1 

•38 


1-8  F8.3  TTF(l) 

9-16  F8.3  TTF (2) 


73*80  F8.3  TTF (10) 


1-8  F8.3 


Continue  same 
as  6a,  6b,  etc. 


TTF(ll) 


ttf(nwl) 


Values  of  transfer  function 
corresponding  to  A/L  values  on 
Cards  6a,  6b,  etc. 


Transfer  function  must  be  in  form 
indicated  by  I FORM  on  Card  3- 


TRANSFER  FUNCTION  - OPTION  3 0 FREQ-3) 

FORMAT  FORTRAN  EXPLANATION 

F8.3  TWMIN  ] Minimum,  maximum,  and  increments  of 

F8.3  TWMAX  Y nondimensionai  encounter  frequency 
F8.3  TWINC  to  .VT77  at  which  transfer  function 

^ is  input 


COLUMNS 

FORMAT 

FORTRAN 

1-8 

F8.3 

TWHIN 

9-16 

F8.3 

TWMAX 

17-24 

F8.3 

TWINC 

1-8 

F8.3 

TTF(l) 

9*16  F8.3  TTF(2) 


Transfer  function  for 


73~80  F8-3  TTF(IO)  I increments  of  u^TfL/g  = TWINC 


1-8  F8-3  TTF(ll)  NTW  « TWHA^;,n^WMIH  + 1 


Continue  with 
10  values  per  card. 
Max.  is  5 cards. 


TTF(NTW)  Transfer  function  at  = TWMAX 


Note:  The  OTNSRDC  Pitch-Heave  Motion  Program  PHM  (formerly  designated 

YF17)  has  a provision  for  punching  pitch,  absolute  motion,  and 
relative  motion  transfer  functions  according  to  the  Option  3 
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TRANSFER  FUNCTION  - OPTION  4 .(IFRECHO 


CARO 

COLUMNS 

FORMAT 

FORTRAN 

EXPLANATION 

6a 

1-8 

F8.3 

YR(1) 

9-16 

• 

9 

F8.3 

• 

• 

YR(2) 

9 

9 

Values  of  wave  frequency  u 
transfer  function  is  input 

9 

73-80 

F8.3 

9 

YR(10) 

l 

Continue  with 
10  values  per  card. 
Max.  is  5 cards 


YR(NWL) 


Must  be  in  ascending  order. 


1-8 

tr\ 

• 

00 

u. 

TTF{1) 

9-16 

F8-3 

TTF(2) 

« 

* 

9 

9 

• 

9 

9 

73-80 

9 

F8.3 

• 

TTF(IG) 

1-8 

F8.3 

TTF(U) 

Continue  with 
10  values  per  card. 
Max.  is  5 cards 


Values  of  transfer  function 
corresponding  to  a values 
on  Cards  6a,  6fc,  etc. 

Transfer  function  must  be  in 
form  indicated  by  1F0RM  cn 
Card  3. 


TTF(NWL) 


4- 


DTNSRDC  ISSUES  THREE  TYPES  OF  REPORTS 

HI  OTNSROC  REPORTS.  A FORMAL  SERIES  PUBLISHING  INFORMATION  OF 
PERMANENT  TECHNICAL  VALUE.  DESIGNATED  BY  A SERIAL  REPORT  NUMBER 

(21  DEPARTMENTAL  REPORTS.  A SEMIFORMAL  SERIES.  RECORDING  informa 
TION  OF  A PRELIMINARY  OR  TEMPORARY  NATURE.  OR  OF  LIMITED  INTEREST  OR 
SIGNIFICANCE.  CARRYING  A DEPARTMENTAL  ALPHANUMERIC  IDENTIFICATION 

(3)  TECHNICAL  MEMORANDA.  AN  INFORMAL  SERIES.  USUALLY  INTERNAL 
WORKING  PAPERS  OR  DIRECT  REPORTS  TO  SPONSORS.  NUMBERED  AS  TM  SERIES 
REPORTS.  NOT  FOR  GENERAL  DISTRIBUTION 
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